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INTRODUCTION 


The  Navigation  Package  (NAVPAC)  is  a  sateilite-bome  data  collection  system  developed  for  the 
Defense  Mapping  Agency  (DMA).  It  consists  of  two  sensor  subsystems: 

1 .  A  receiver/data  subsystem  capable  of  tracking  up  to  three  Navy  Navigation  Satellites 
(NAVSATs)  simultaneously  and  recording  integrated  Doppler,  first-order  ionospheric  refraction, 
and  timing  information. 

2.  A  miniature  electrostatic  accelerometer  (MESA)  capable  of  measuring  accelerations  due  to 
atmospheric  drag,  winds,  and  orbit  adjust  thrusts. 

The  receiver/data  subsystem,  which  also  provides  the  power,  command,  telemetry,  and  MESA  data 
storage  functions  was  designed  and  built  by  the  Applied  Physics  Laboratory  of  Johns  Hopkins 
University  (APL/JHU)  according  to  the  specifications  provided  by  the  Naval  Surface  Weapons 
Center  (NSWC).  The  MESA  was  designed  and  built  by  Bell  Aerospace  TEXTRON  according  to  the 
specifications  provided  by  the  Air  Force  Geophysics  Laboratory  (AFGL)  and  NSWC. 

The  NAVPAC  satellite-to-satellite  tracking  (SST)  Doppler  data  and  the  accelerometer  data  are 
used  by  NSWC  for  orbit  determination  for  the  host  satellite.  In  the  future,  these  data  may  also  be 
used  for  gravity  field  model  improvement.  The  MESA  data  are  also  used  by  AFGL  for  atmospheric 
density  and  wind  studies. 

NAVPAC  stores  the  collected  data  in  a  core  storage  unit  until  it  becames  full  (approximately 
one  rev).  At  that  time,  the  data  are  automatically  transferred  to  an  onboard  tape  recorder.  Daily, 
the  data  accumulated  on  this  tape  recorder  are  telemetered  to  a  ground  station,  recorded  on  analog 
tape,  and  provided  to  NSWC.  Here,  the  analog  tape  is  digitized,  preprocessed  in  the  “Quick  Look” 
mode  using  preliminary  trajectories,  and  a  preliminary  data  quality  evaluation  is  done.  Later,  the 
data  are  grouped  for  a  given  time  span  and  preprocessed  in  the  “Normal”  mode  using  the  precise 
NAVSAT  trajectories.  After  preprocessing,  the  data  are  combined  with  the  station  tracking  data 
for  the  host  satellite  in  NSWC’s  orbit  computation  program.  CELEST,  to  produce  the  required 
epheme  rides. 

The  purpose  of  this  report  is  to  provide  a  detailed  mathematical  description  of  the  NAVPAC 
satellite-to-satellite  tracking  data  processing  procedures  from  preprocessing  to  orbit  determination 
including  an  overview  of  the  hardware  operation  required  for  the  correct  interpretation  of  the 
Doppler,  refraction,  and  timing  data.  A  similar  description  for  the  NAVPAC  accelerometer  data  is 
given  in  Reference  I . 


RECEIVER/DATA  SUBSYSTEM  OVERVIEW 


JJie  NAVPAC  receiver/data  subsystem  includes  three  identical,  independent  receivers  for 
phase-locked  loop  tracking  of  the  coherent  400/1 50-MHz  signal  pair  broadcast  by  the  Navy  Naviga¬ 
tion  Satellites.  Each  receiver  is  similar  in  design  to  the  ground-based  Geoceiver  tracking  equipment. 
Both  are  designed  in  the  “master/slavc”  configuration:  i.e.,  the  scaled  signal  dynamics  from  the 


400-MHz  channel  is  subtracted  from  the  150-MHz  channel  leaving  only  signal  variations  due  to  the 
difference  in  ionospheric  refraction  effects  on  the  two  signals.  The  low  channel  can  track  a  satellite 
only  when  the  high  channel  is  phase-locked. 

The  signal  dynamics  that  the  NAVPAC  receivers  must  track  are  over  twice  what  a  ground- 
based  Geoceiver  would  see.  The  maximum  Doppler  shift  is  ±  20  kHz,  the  maximum  Doppler  rate  is 
400  Hz/sec,  and  the  maximum  Doppler  acceleration  is  ±  6  Hz/sec2.  Figure  1  shows  worst-case  plots 
for  these  three  parameters  for  a  hypothetical  satellite-to-satellite  pass.  Since  the  actual  NA  VSAT 
broadcast  frequencies  are  offset  by  -80  ppm  (-32  kHz)  at  400  MHz,  but  the  local  NAVPAC  oscil¬ 
lator  is  not  offset,  the  resulting  Doppler  frequencies  range  from  1 2  to  52  kHz.  Acquisition  always 
occurs  in  the  1 2-  to  32-kHz  range  at  elevation  angles  no  greater  than  -  10°  as  viewed  from  the  host 
satellite.  Each  receiver  sweeps  in  this  frequency  range  asynchronously  with  the  others  with  each 
sweep  taking  40  sec  (20  kHz  divided  by  500  Hz/sec).  If  a  signal  is  detected  by  the  high  channel, 
the  sweep  stops  and  NAVSAT-type  modulation  is  searched  for.  If  it  is  not  found,  the  sweep  is 
restarted  and  a  bogie  counter  is  incremented.  If  the  correct  modulation  is  present,  an  indicator  is 
sent  to  the  data  subsystem  to  initiate  the  Doppler  count  and  to  the  low  channel  to  allow  it  to  begin 
tracking. 

The  low  channel  actually  tracks  the  frequency  given  by 
fL  =  4(eL-^eH> 

where 


ei  =  error  in  frequency  at  1 50  MHz  due  to  refraction 
eji  =  error  in  frequency  at  400  MHz  due  to  refraction 

To  first  order,  the  error  due  to  refraction  is  inversely  proportional  to  the  frequency.  Therefore 
eL  =  8/3  e^-  Substituting  this  into  the  above  equation  gives 

3  „  8  3  1 

fL  =  4(eL  -  g  eH  >  =  4(^eH  -  geH  >  =  9^H 

Thus,  the  low-channel  tracked  frequency  is  approximately  9-1/6  times  the  error  in  frequency  at 
400  MHz  due  to  refraction.  The  low  channel  can  track  frequencies  less  than  ±  600  Hz  with  a  maxi¬ 
mum  rate  of  50  Hz/sec.  The  refraction  count  takes  place  in  the  data  subsystem.  If  a  loop  open  and 
sweeping  indicator  is  received  by  this  subsystem  at  any  time  during  the  count,  the  count  is  stopped 
and  the  refraction  count  for  this  measurement  is  set  to  zero.  (In  the  original  NAVPAC  design  the 
loop  not  phase-locked  indicator  was  used  instead.) 

Extensive  priority  logic  is  built  into  the  receiver  subsystem  controller.  If  one  receiver  is  locked 
onto  and  tracking  a  N  A  VS  AT  signal  and  another  receiver  encounters  the  same  signal  during  its 
normal  sweep  mode,  the  sweep  stops  and  the  Doppler  frequencies  are  compared.  If  they  are  the 
same  (within  25  Hz),  the  original  receiver  has  priority  and  continues  to  track  while  the  other 
receiver  resumes  its  sweep.  Whenever  phase-lock  is  lost  for  6  sec,  the  receiver  reverts  to  the  sweep 
mode  to  try  to  reacquire  the  signal.  It  now  sweeps  over  a  frequency  range  starting  at  fy  -  1 .4  kHz 
(where  fy  is  the  frequency  at  which  loss  of  lock  occurred)  and  extends  20  kHz  or  to  52  kHz  which¬ 
ever  is  less.  This  sweep  range  is  chosen  because  in  certain  cases  positive  Doppler  rates  are  possible. 


Wo. 


Figure  2  shows  the  dynamics  of  a  hypothetical  satellite-to-satellite  pass  with  large  positive  Doppler 
rates.  This  sweep  continues  until  the  signal  is  reacquired  or  four  times  at  which  time  the  receiver 
resets  its  sweep  range  back  to  12  to  32  kHz.  It  is  possible  that  one  of  the  other  two  receivers 
could  lock  on  to  this  signal  before  the  original  receiver  reacquires  it.  In  this  case  the  new  receiver 
has  priority.  Another  important  case  is  when  two  receivers  are  tracking  separate  signals  and  their 
Doppler  curves  cross.  Several  results  are  possible.  The  receivers  may  continue  to  track  their  original 
signals.  One  receiver  may  switch  to  the  other’s  signal  in  which  case  the  receiver  with  priority  (de¬ 
termined  by  the  order  in  which  they  started  tracking)  continues  to  track  the  common  signal  and 
the  other  reverts  to  the  reacquisition  sweep  mode.  It  may  reacquire  the  remaining  signal  present  or 
the  third  receiver  may  acquire  this  signal.  Another  possible  outcome  is  that  the  receivers  will  just 
switch  signals  without  any  loss  of  lock.  In  all  cases  the  priority  logic  circuitry  ensures  that  the 
same  signal  is  never  tracked  by  more  than  one  receiver. 

This  logic  can  be  disabled  by  command  allowing  all  three  receivers  to  track  the  same  satellite. 
This  is  called  the  priority  defeat  mode.  In  this  mode  the  Doppler  count  intervals  are  not  simul¬ 
taneous.  This  is  because  the  sweeps  are  asynchronous  and  acquisition  is  independent  for  each  re¬ 
ceiver.  Commands  also  exist  for  turning  off  individual  receivers  and  disabling  the  modulation  search 
circuitry. 

Timing  for  the  Doppler  data  is  based  on  a  master  clock  frequency  of  833.333  kHz  derived 
from  an  ultrastable  5-MHz  quartz  crystal  oscillator.  This  clock  frequency  corresponds  to  a  time 
resolution  of  1 .2  psec.  The  time  code  generator  (TCG)  is  then  just  a  42-bit  counter  that  increments 
every  1 .2  psec  and  starts  when  the  data  subsystem  is  commanded  on.  This  counter  recycles  to  0 
approximately  every  6 1  days  or  when  reset  by  command.  Whenever  the  data  subsystem  receives 
a  high-channel,  phase-locked  indicator  from  a  receiver,  the  Doppler  counter  starts  at  the  first 
negative  zero  crossing.  Cycles  are  counted  for  exactly  30  NAVPAC  sec  (25  X  10' 2  TCG  counts). 

At  this  time  a  time  over  counter  (TOC)  starts  counting  at  1 .2  psec  steps  until  the  next  negative 
zero  crossing  occurs.  Also  at  this  time  the  TCG  is  read  out  and  put  in  the  data  stream  along  with 
a  receiver  identification  (ID)  and  is  identified  as  a  Doppler  time  mark  (DTM).  (The  TCG  is  not  read 
out  at  the  beginning  of  the  count  just  in  case  the  count  is  not  completed.)  The  Doppler  count 
register  now  contains  one  less  than  the  total  number  of  Doppler  counts  collected  in  an  interval  of 
30  +  TOC  X  1 .2  X  10~6  NAVPAC  seconds.  A  Doppler  word  (DW)  is  then  placed  in  the  data  stream 
containing  the  Doppler  count,  time  over  count,  and  receiver  ID.  At  the  end  of  the  time  over  count, 
the  Doppler  counter  restarts  at  zero  and  this  sequence  is  continued  as  long  as  phase-lock  is  main¬ 
tained.  Doppler  time  marks  are  placed  in  the  data  stream  every  8  Doppler  words  as  long  as  the 
receiver  has  not  lost  lock. 

In  parallel  with  the  Doppler  count  is  a  refraction  count  that  terminates  at  the  end  of  the 
30-sec  portion  of  the  Doppler  count.  (An  error  of  1  refraction  cycle  corresponds  to  an  error  in 
the  range  difference  measurement  of  approximately  8  cm.)  The  Doppler  word  mentioned  above 
also  contains  the  refraction  count. 

All  of  the  Doppler  timing  is  an  internal  NAVPAC  function  based  on  the  output  of  the  local 
oscillator  driving  the  TCG.  The  exact  start  time  of  the  master  clock  and  the  exact  frequency  at 
any  instant  are  not  known.  To  allow  for  conversion  from  NAVPAC  time  to  absolute  time  (UTC) 
in  software.  NAVPAC  contains  a  navigation  message  recovery  (NMR)  unit.  This  unit  demodulates 
the  signal  to  get  the  navigation  message  on  the  400-MHz  channel  and  decodes  the  two-minute  time 
mark  (TMTM)  and  satellite  ID.  These  NAVSAT  TMTMs  are  kept  synchronized  with  UTC  very 
accurately.  The  NMR  unit  cycles  among  the  receivers  until  a  Doppler  lock  indication  is  present. 


It  then  decodes  the  signal  in  order  to  recover  the  TMTM  and  satellite  ID.  If  4  min  have  elapsed 
without  a  successful  recovery,  the  NMR  cycles  to  the  next  locked  receiver.  When  3  TMTM  is 
decoded,  the  least  significant  3b  hits  of  the  TC  Ci  register  along  with  the  satellite  and  receiver  IDs 
are  placed  in  the  data  stream  identified  as  a  TMTM.  This  portion  of  the  TCG  register  recycles  to  0 
approximately  every  2d  hr. 

TMTMs  are  the  only  link  between  a  receiver  and  the  satellite  it  is  tracking.  It  is  possible  for  a 
TMTM  to  occur  in  the  data  stream  before  a  DTM  signifying  acquisition  if  the  TMTM  occurs  during 
the  first  30-sec  Doppler  count.  A  TMTM  may  also  occur  during  a  Doppler  count  that  is  not  com¬ 
pleted  due  to  loss  of  lock.  It  is  also  possible  for  a  TMTM  from  a  particular  satellite  to  be  obtained 
by  one  receiver  and  actual  acquisition  of  this  satellite  (defined  by  the  first  complete  Doppler 
count)  may  be  on  another  receiver.  If  three  satellites  are  being  tracked  simultaneously,  as  much  as 
6  min  can  elapse  between  indications  that  each  receiver  is  still  tracking  the  same  satellite. 

DTMs,  DWs.  and  TMTMs  are  placed  in  the  data  stream  as  they  are  collected  regardless  of  which 
receiver  they  come  from  or  which  satellites  are  being  tracked.  NAVPAC  also  has  the  capability  to 
time  tag  three  external  events.  All  of  the  corresponding  data  words  are  48  bits  as  described  in 
Figure  3. 
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Figure  3.  Format  of  the  4H-Bit  Data  Word 


The  NAVPAC  data  subsystem  also  collects  both  analog  and  discrete  telemetry  from  the  re¬ 
ceivers.  The  analog  data,  which  are  digitized  into  X-bit  words,  consist  of  AGC  readings  for  all 
receivers  (high  and  low  channels)  and  are  sampled  twice  each  telemetry  frame  span  of  1  14,5  see.  The 
discrete  data  consist  of  such  quantities  as  Doppler  lock  and  sweep  status  indicators  for  all  receivers 
(high  and  low  channels)  and  bogie  counters  (breach  high  channel.  These  are  sampled  once  every 
1  14.5  sec.  NAVPAC  also  lias  a  troubleshooting  telemetry  collection  mode  in  which  data  are  sampled 
at  eight  times  the  normal  rate. 

The  above  discussion  was  only  intended  to  give  the  reader  a  sufficient  understanding  of  the 
receiver  data  subsystem  operation  for  correctly  interpreting  the  NAVPAC  Doppler,  refraction, 
and  timing  data.  Reference  2  contains  a  complete  detailed  description  of  the  NAVPAC  hardware 
built  by  APL/JHU. 


DATA  PREPROCESSING  (NAVPAC) 


The  primary  purpose  of  preprocessing  the  raw  NAVPAC  Doppler  data  (using  the  NAVPAC 
computer  program)  is  to  convert  it  into  a  form  usable  by  the  QSOLVK  data  quality  evaluation 
program  and  by  the  CEL. I  ST  orbit  computation  program.  In  addition,  information  about  the 
NAVPAC  clock  is  generated  for  use  in  correctly  time-tagging  the  telemetry  and  accelerometer  data. 
Preprocessing  involves  four  passes  through  the  data  to  accomplish  the  following: 

1 .  Calibrating  the  clock 

2.  Sorting  data  into  passes 

3.  Time-tagging  and  converting  counts  to  range  differences 

4.  Reordering  passes  by  NAVSAT 

The  raw  data  are  input  as  a  stream  of  48-bit  words  as  defined  in  the  previous  section.  I  wo  modes  of 
preprocessing  are  possible:  “Quick  Look”  and  “Normal.”  These  are  identical  except  that  the  data 
for  only  a  subset  of  the  possible  NAVSATs  are  processed  in  the  “Quick  Look”  mode. 


CALIBRATING  THE  CLOCK 

As  mentioned  in  the  receiver/data  subsystem  overview,  the  TMTMs  relate  NAVPAC  time 
(contents  of  a  counter  incrementing  every  1 .2  NAVPAC  /usee)  to  Coordinated  Universal  Time 
( UTC).  The  input  NAVPAC  clock  reset  time  is  only  an  approximate  UTC  time  at  which  the  data  sub¬ 
system  on  command  was  executed.  This  uncertainty  along  with  the  fact  that  the  time  increment 
is  not  exactly  1 .2  UTC  /usee  and  varies  continuously  imply  that  some  type  of  adjustment  needs  to 
be  made  to  (lie  NAVPAC  times  ( DTMs,  host  vehicle  time  mrrks  (HVTMs),  and  telemetry  and 
accelerometer  time  marks)  based  on  the  received  1  MTMs. 

To  ensure  that  the  displayed  times  for  all  the  printed  raw  Doppler  data  can  later  be  compared 
to  calibrated  times  for  analysis  purposes, the  reset  time  is  recomputed  based  on  the  first  TMTM  for 


which  trajectory  information  is  available.  The  following  definitions  are  based  on  the  convention 
that  r  represents  a  relative  time  and  t  an  absolute  time: 


’TMTM 

AtTT 

A,l)l> 

A,RI> 


NAVPAC  time  of  TMTM  in  counts  X  1 .2  X  10"6  see 

*  |rNS  ~  rllS | 

transmission  time  = - =  - — - 

c  c 

decode  delay 
ith  receiver  delay 

input  reset  time 


(K  +  tTMTM  ’s  ,usl  founded  to  the  nearest  even  minute  mark.  fXMTM •  The  approximate  NAVPAt 
time  at  the  time  of  emission  of  the  I  MTM  by  the  satellite,  .  is  then  given  by 


tTMTM  -  rTMTM  '  AtTT  Atl)i>  "  AtRI)j 
1  he  recomputed  reset  time  is  then  given  by 
tRrei,.„,pu.«l  =  ‘™™  "rTMTM 

In  addition,  the  next  two  I  M  l  Ms  are  compared  against  the  first  to  verity  that  the  first  was  not 
grossly  in  error.  This  recomputed  reset  time  is  used  in  computing  all  times  displayed  in  the  raw  data 
printout. 

In  the  first  pass  through  the  data,  quantities  for  each  TMTM  needed  for  the  clock  calibration 
are  computed  as  given  later  in  this  section.  Also  if  the  NAVPAC  clock  recycles  or  has  recycled, 
the  DTMs  and  UVTMs  are  adjusted  accordingly  where  the  increment  is  integer  multiples  of 
24-  =  439X04051  1 104  (approximately  bl  days  2  hr).  Since  only  36  bits  of  the  42-bit  time  word  are 
present  in  the  TMTMs.  each  I'M  I'M  count  is  compared  against  the  previous  DTM  and  adjusted  in 
increments  of  2-,<1  =  6X719476730  (approximately  23  hr)  until  larger.  All  TMTMs.  DTMs.  and 
UVTMs  are  given  an  associated  day,  hour,  and  minute  label  and  every  data  word  (including  Doppler 
words)  is  given  a  time  tag  in  seconds  from  the  beginning  of  the  day  in  the  raw  data  printout.  The 
D  I  M  times  printed  have  30  sec  substraeted  from  them  to  indicate  the  start  time  of  the  following 
Doppler  count.  I  or  this  reason  TMTMs  and  DTMs  may  appear  to  be  out  of  order  in  the  printout. 
The  time  tag  associated  with  each  l)W  is  the  end  time  of  that  particular  count  to  the  nearest 
second.  The  begin  time  of  a  particular  count  is  either  the  end  time  of  the  last  count  (on  the 
same  receiver)  or  the  time  of  the  preceding  DTM.  The  refraction  count  present  in  each  DW  has  had 
2IS  =  3276X  subtracted  from  its  raw  value  to  get  the  true  refraction  count  (which  may  be 
positive  or  negative ). 

To  relate  the  satellite  ID  in  the  TMTM  data  word  to  the  NSWC  satellite  numbers  a  table 
look-up  procedure  is  used,  fable  I  gives  the  correspondence  between  TMTM  IDs.  NSWC  satellite 
numbers,  and  API  satellite  numbers  for  all  current  NAVSATs.  Hie  printed  IDs  are  the  NSWC 
satellite  numbers  unless  no  match  is  found  in  which  case  the  original  recovered  ID  is  printed  and 


X 


the  TMTM  is  deleted  from  any  further  processing.  The  total  number  of  bad  IDs  is  printed  after  the 
pass  summary  table  to  be  described  later. 

fable  1.  Satellite  Numbers 


TMTM  ID  NSWC  APL 


36 

58 

30120 

40 

59 

30130 

56 

60 

30140 

28 

68 

30190 

16 

77 

30200 

52 

64 

30180 

62 

86 

30460  (TIP  11) 

63 

90 

30470  (TIP  III) 

32 

93 

301 10  (TRANSAT) 

After  identifying  the  satellite  from  which  the  TMTM  was  received,  tR  ted  +  tTmtM  is 
rounded  to  the  nearest  even  minute  mark,  t-j-j^TM  •  Next  tlle  tilne  °*  transmission,  Atpp,  is 
computed  as 


p(tTMTM>  =  lrNS~  rHsl 
c 

NAVSAT  position  at  tp^pM  /  *TMTM  's  assumed  to  be  a  trajectory 
host  satellite  position  at  tp\jpM  (  timeline  for  each  satellite 

Then  the  approximate  UTC  time,  tj^j-pM  •  at  which  the  TMTM  was  received  by  NAVPAC  is  given  by 

lTMTM  =  lTMTM  +  AtTT  +  A,DI)  +  AtRD. 

1 

where  Atpp  =  decode  delay  =  .007373  sec  (the  same  for  each  receiver) 

and  AtRp  =  i1^1  receiver  delay  (different  for  each  receiver)  i  =  1 .  2,  or  3 

This  is  approximate  because  none  of  the  three  terms  added  to  t-p^TM  is  known  exactly,  e.g..  AtRp> 
is  actually  a  function  of  signal  strength.  The  above  computation  is  repeated  for  each  TMTM  foi  1 
which  Atpp  can  be  computed  resulting  in  two  values  for  each  time  mark  rpMTM  and  tp^p^ .  This 
correspondence  between  NAVPAC  and  UTC  times  is  also  in  error  because  the  NAVSAT  clocks  are 
not  synchronized  exactly  with  UTC  or  each  other. 

A  polynomial  equation  is  used  to  convert  from  a  NAVPAC  time,  r,  to  a  UTC  time,  t,  as 
follows: 


TT 


where  r^s  = 
rHS  = 


t  =  tR  +  A  +  (l  +B)t+  Ct2 

^recomputed 


(i » 


9 


This  can  be  rewritten  as: 


recomputed 


-  t  =  A  +  Br  +  Ct* 


For  each  TMTM,  an  approximate  value  for  t  is  computed  giving  an  equation  of  the  form  : 


*TMTM.  “  lR 


recomputed 


tTMTM;  “  A  +  BtTMTM.  +  (  tTMTM; 


j  =  I.  2,  ...  N 


TMTM 


This  is  an  overdetermined  system  of  linear  equations  in  the  unknowns  A,  B,  and  C  and  can  be 
solved  using  least-squares  techniques  without  a  priori  information.  The  normal  equations  are  given 
by: 


TMTM 


TrTMTM;  ?rTMTM 


JrTMTMj  yTTMTM)  ^rTMTMj 
frTMTM  ?rTMTM.  ?tTMTM, 


lTMTM  '  lTMTM 


k  recomputed 


'TMTM. 


Vt" 

j  mMTMj 

}ttTMTMj  ‘TMTMj 
?rTMTMj  tTMTMj 


and  all  sums  are  from  j  =  1  to  NTM 

After  these  eq1  ations  are  solved  for  A,  B.  and  (\  the  individual  residuals  of  fit,  €j,  are  computed 
where 


cj  =  (tTMTM  "  *R 


recomputed  “  rTMTM/-  (A  +  BrTMTMj  +  CtT MTM}> 


The  root  mean  square  ( RMS)  of  these  residuals  is  then  computed 


'TMTM 


and  all  TMTM  for  which  |ej|  >  nto[  crMS  are  deleted  and  the  least -squares  solution  repeated. 

The  above  iterative  procedure  is  done  until  no  more  TMTMs  are  deleted  or  until  the  maximum 
number  of  iterations  is  reached.  For  each  iteration,  the  NAVPAC  time  count  and  the  correspond¬ 
ing  residual  for  all  deleted  points  are  printed.  For  the  last  iteration,  the  calibrated  TMTMs,  a 
correction  term,  the  t’rVITM j  terms,  and  the  final  residuals  are  printed  along  with  the  original  and 
final  number  of  TMTMs.  the  number  of  iterations,  the  RMS  of  the  residuals  (in  p sec),  the  re¬ 
computed  reset  time,  and  the  A.  B,  and  l  values.  The  correction  term  mentioned  above  is  given 
by  BfTMTMj  +  O^MTM,  This  is  the  clock  correction  due  only  to  the  oscillator  frequency  offset 
(related  to  B)  and  frequency  drift  (related  to  C). 

liquation  ( I )  with  tR„  ,  (ed.  A,  B,  and  C  as  determined  above  is  then  used  to  convert  any 
NAVPAC  time,  r.  to  its  equivalent  UTC  time,  t.  After  the  clock  calibration  but  still  in  this  first 


10 


pass  through  the  d->»a.  the  lIVTMs.  which  have  already  been  separated  hy  type,  are  calibrated  and 
the  appropriate  delays  (different  tor  each  type!  subtracted  out 

SORTING  DATA  INTO  PASSE  S 

Because  of  simultaneous  tracking  of  up  to  three  NAVSATs.  the  possible  receiver  switching 
that  can  occur,  anil  the  loss  of  lock  and  reacquisition  possibilities,  sorting  of  the  Doppler  data  into 
passes  is  a  complicated  procedure.  I  he  key  data  words  in  the  sorting  are  the  I  M  I  Ms.  since  they 
provide  the  only  connection  between  receiver  number  and  satellite  II)  All  data  for  legitimate  satel¬ 
lite  IDs  are  sorted.  I  he  following  description  gives  a  set  of  rules  for  assigning  data  to  passes  based  on 
DTMs  and  TMTMs. 

If  a  DTM  occurs  more  than  10  min  after  the  last  DEM  from  any  receiver,  it  is  assumed  that  all 
previously  active  passes  have  ended  am  I  a  possible  new  pass  has  begun,  lour  DWs  are  required 
before  a  new  pass  is  declared  active.  If  a  DIM  occurs  on  a  receiver  that  has  been  inactive,  the 
assumption  is  made  that  this  indicates  the  start  of  a  new  pass.  If  the  DTM  is  the  second  DTM  from 
a  receiver  without  an  intervening  I'MTM.  a  pass  has  started  but  the  satellite  has  not  been  identified. 
If  a  DTM  and  a  TMTM  have  already  been  found  on  this  receiver,  there  are  two  cases  to  consider: 

I )  If  the  time  since  the  last  DTM  is  exactly  4  NAVPAC  minutes  +  TOCj.  all  DWs  between  these 
D  I  Ms  are  assigned  to  the  current  pass.  2)  If  this  time  difference  is  not  exact,  future  TMTMs  arc 
checked  to  see  if  a)  the  receiver  lost  the  signal  but  reacquired  it  later  in  which  case  the  pass  is 
continued  or  b)  a  different  satellite  ID  is  found  in  which  case  the  DTM  defines  the  beginning  of  a 
new  pass.  Also  in  case  bl  the  previous  pass  is  not  terminated  because  the  possibility  still  exists  that 
this  satellite  may  be  reacquired  by  another  receiver. 

When  a  TMTM  is  being  processed,  its  receiver  number  and  satellite  ID  are  compared  against 
previous  TMTMs  until  at  least  a  partial  match  is  found  or  ID  min  have  elapsed,  four  outcomes  are 
possible: 

1 .  Receiver  and  satellite  match  f his  indicates  that  between  the  two  TMTMs  in  question  all 
the  data  from  this  receiver  are  from  the  same  satellite.  This  is  normally  the  case. 

2.  Receiver  and  satellite  do  not  match  I  his  gives  no  information  on  the  present  receiver  so 
the  next  I  M  l  M  back  is  checked.  If  no  match  is  found,  it  is  assumed  that  the  current  TMTM  indi¬ 
cates  the  beginning  of  a  new  pass 

.V  Receiver  matches,  satellites  do  not  This  indicates  that  the  receiver  was  previously  track¬ 
ing  another  satellite.  I  he  present  TM  TM  then  may  indicate  the  beginning  of  a  new  pass  or  a  receiver 
switch  The  satellite  that  was  previously  tracked  is  searched  for  later  on  another  receiver. 

4.  Satellite  matches,  receivers  do  not  11ns  also  indicates  a  switch  The  pass  from  the  pre¬ 
vious  receiver  must  be  kept  active  in  case  its  satellite  ID  reappears  later. 

When  a  satellite  switch  occurs  without  the  involved  receiver  losing  lock,  an  ambiguity  exists  as  to 
during  exactly  which  Doppler  count  the  switch  took  place.  It  is  assumed  that  all  DWs  since  the  last 
TMTM  come  from  the  most  recent  satellite  tracked  by  this  receiver. 


l-'or  data  collected  in  the  priority  defeat  inode,  the  sorting  is  simplified.  All  tests  for  switch¬ 
ing  are  deleted,  each  receiver  is  treated  independently,  and  the  only  requirement  is  to  determine 
when  a  particular  pass  ends  and  the  next  one  begins. 

The  main  result  of  this  pass  through  the  raw  SSI  Doppler  data  is  the  assignment  of  pass 
numbers  and  the  corresponding  satellite  ID  to  every  DI  M  and  D\V 


TIME  TAGGING  AND  CONVERTING  COUNTS  TO  RANGE  DIFFERENCES 

The  data  are  now  reprocessed  in  order  to  assign  calibrated  times  to  each  observation  and  to 
convert  from  uncorrected  Doppler  counts  at  4HI)  Mil/  and  rvliuction  counts  to  I  list  order  refraction 
corrected  range  difference  data  with  a  nominal  frequency  luas  removed.  In  addition,  retraction 
counts  are  examined  for  zero  or  large  values  in  which  case  the  data  are  tagged. 

I  ach  pass  begins  with  a  DIM.  Since  the  time  tag  associated  with  a  range  difference  observa¬ 
tion  is  the  end  time  of  the  corresponding  Doppler  count,  a  dummy  observation  must  he  included  at 
the  beginning  of  each  pass  which  specifies  the  start  time  of  the  first  count  I  his  dummy  observation 
consists  of  setting  the  observed  v.,uie  to  O.atnl  the  observation  sigma  to  ItHHL  1  he  I  K  time,  t. 
corresponding  to  a  D  I  M  is  given  by 


‘DIM  =  ,rf)TM  30)eal 


I  his  notation  (  )CJ|  means  substitution  of  the  NAVPAC 
time  within  (  >  into  the  time  calibration  liquation  (1 ). 


I  lie  ETC  time  associated  with  the  first  observation  immediately  following  the  DI  M  is  given  by 


'obs.  =  <TI)  I  M  f  IO(  1  'cal 


Therefore,  the  UK  time  associated  with  the  i1'1  observation  trom  this  same  receiver  given  the 
NAVPAC  time  associated  with  the  i  -  Is*  observation  is  given  by 


'uhs 


>bs  ,  +  *>  +  <<><i>cal 


(2) 


where 


ohs 


r  I )  1  M 


+  li  I 


l  ■  40  a  roc,  +  ...  a  |<  >C 


In  all  these  expressions,  the  lot  is  actually  the  IOC  front  the  DW  multiplied  In  12  x  10  h  sec. 

I  lus  procedure  continues  until  the  next  DIM  is  encountered  for  the  current  pass.  If  it  is  from  the 
same  receiver.  X  I )V\ ->  Irom  this  receiver  have  been  processed  since  the  last  DIM.  and 


rl)TM. 


rl)l  M, 


+  240  a  1  IOC 

i  1 


then  no  loss  ol  loch  has  occurred  and  theretore  liming  is  continued  based  on  I  quation  ( 2 ).  In  all 
other  cases,  timing  is  restarted  based  on  this  new  I)  I  M  with  a  dummy  observation,  l  ime  tagging  for 
up  to  three  passes  can  he  taking  place  concurrently  .  I  lie  only  timing  ambiguity  exists  when  a  pass 
switches  from  one  receiver  to  another  without  a  loss  of  loch  (no  DI  M).  In  this  case,  the  observation 
times  must  he  based  on  the  next  available  DIM  from  this  receiver  Ivy  subtracting  40  +  IOC,  from 
each  uncalibrated  observation  lime  starting  with  the  first  observation  alter  t he  DIM  and  then 


calibrating  the  recall  All  l'  l(  observation  times  are  converted  to  a  day  number  am!  seconds  from 
the  beginning  of  this  day  alter  calibration 


I  or  each  I  trappier  word,  the  Doppler  count  « 1 M  ,  I  and  re  tract  ion  count  ( RC  j )  are  com¬ 
bined  as  follows  to  get  the  range  difference  observation 

e  f  0 

Kl),  =  -  -  .  (D(  ,  ..  K<  >  f  ,  t  30  +  |()(  » 

'  '  +  t  Offset  L  '  » 

where  DC  =  Doppler  count  at  4(10  Mil/  (  raw  Doppler  count  +  1  if'KX'jX).) 

RC  =  refraction  count  (raw  retraction  count  -  32768) 

I  _  \  d2()0()  11/ i  80  ppm  at  400  Mil/)  for  satellites  58,  5{).  60,  64. 68,  and  77 

°^'scl  j  33702(1/1  84.48  ppm  at  400  Ml l/t  for  all  other  satellites 

30.+  fO(  ’  =  Doppler  count  interval  in  NAVPAf  seconds 

c  =  speeil  of  light  (in  km  sect 

f  =  nominal  N.M  I’At  frequency  =  400  X  I0(’  Hz 

I  ach  observation  is  assigned  an  observation  sigma  of  ( 30.  +  TOCjf  1  “  ■  However,  if  R(  j  =  0  or 
RC  ,!  •  2000.  the  sigma  is  set  to  1000..  Also,  it  i R Dj  >  500.  km.  the  observation  value  is  set  to  0. . 

After  all  the  D  I  Ms  and  DVVx  for  a  given  pass  are  processed,  the  pass  is  assigned  a  type  number 
based  on  the  angular  separation  of  the  host  satellite  and  NAVS.A  f  orbit  planes.  Table  2  defines  the 
type  numbers  as  a  function  of  the  approximate  right  ascension  of  the  ascending  node  difference. 

Pass  ty  pes  2  and  3  are  usually  short  (less  than  I  3  of  a  rev  i  and  have  a  unique  time  of  closest 
approach  I  1(  A)  with  range  differences  negative  before  and  positive  after  this  time.  An  approximate 
l<  A  is  defined  as  the  firs*  observation  time  for  which  the  observation  changes  sign.  If  only  a  portion 
of  the  actual  satellite  pass  is  present  with  all  observations  of  the  same  sign,  then  the  approximate 
If  A  is  defined  as  the  last  observation  time,  if  all  observations  are  negative  and  the  first  observation 
time,  if  all  observations  are  positive 


I  able  2 


Type  Number 


Magnitude  of  Right  Ascension  Difference  <Ci  I 


I 

s 

3 


0  to  45 
45  to  135 
135  to  181) 


Type  1  passes  are  much  longer  than  the  other  two  types  (up  to  2  hr)  and  may  have  more  than 
one  real  H  A.  I  bese  passes  are  usually  segmented  and  each  segment  treated  as  a  separate  pass  in 
further  processing.  An  example  of  the  Doppler  curve  for  a  hypothetical  pass  of  this  type  was  given 
in  f  igure  2  A  plot  of  the  range  difference  values,  which  are  actually  used  in  the  segmenting 


procedure,  would  have  a  similar  shape  but  with  the  signs  reversed.  Note  that  there  are  long  periods 
when  the  Doppler  shift  is  increasing  (range  differences  are  decreasing).  file  segmenting  starts  by 
searching  the  range  difference  values  for  a  change  insign  or  35  min  whichever  comes  first.  A  change 
in  sign  from  to  +  is  a  normal  I  (  A.  A  change  from  +  to  is  also  called  a  TCA  hut  is  actually  a  time 
of  furthest  approach  (zero  Doppler  shift  hut  a  local  maximum  instead  a  local  minimum  range).  If 
no  TCA  is  found,  the  pass  is  segmented  at  this  time  and  a  type  number  6  assigned  to  the  pass  along 
with  a  1C  A  defined  hy  -  (ti^,  ,)hs  -  tt'irit  0bs. >  -  ' ,lc  negative  sign  attached  to  this  value  indicates 
that  it  is  not  a  true  TCA.  If  a  TCA  is  found,  the  next  local  minimum  or  maximum  or  a  time  span 
of  35  min  defines  the  end  of  this  segment.  If  the  1CA  corresponds  to  a  to  +  sign  change,  type 
number  4  is  assigned  to  the  pass.  If  it  corresponds  to  a  +  to  sign  change,  type  number  5  is  assigned. 
If  a  pass  segment  as  defined  above  is  less  than  10  min  long,  the  segment  is  extended  to  the  next 
local  extremum  or  another  25  min.  If  the  last  observation  in  the  pass  is  reached  and  the  current 
segment  is  less  than  10  min  long,  this  segment  is  combined  with  the  previous  segment.  Finally,  if 
the  last  observation  is  reached  without  any  segmenting  being  done,  the  pass  type  is  left  at  1  and 
again  -  (t^,  0hs  -  t  first  uhs.  *  -  ’s  assigned  as  the  TCA.  All  of  the  above  results  in  pass  segments  less 
than  45  min  in  duration  If  the  segmenting  takes  place  during  an  observation  span  where  no  loss  of 
lock  occurs,  the  first  observation  of  the  new  segment  must  lie  a  dummy  observation  with  the 
observation  time  set  equal  to  that  of  the  last  observation  in  the  previous  segment. 

Various  quantities  designed  to  provide  pertinent  information  about  each  NAVPAC  pass  are 
computed  anti  entered  in  a  pass  summary  table.  A  subset  of  these  quantities  are  also  totalled  over  all 
passes.  Hlevation  and  azimuth  angles  relative  to  the  host  satellite  at  rise.  TCA.  and  set  are  computed 
for  each  pass  after  a  more  exact  TCA  has  been  determined.  TCA  is  refined  by  iterating  the  follow¬ 
ing  formula 

TCAj+j  =  TCAj  ~  (p1  pVK  Aji/ip1  p  +  pF  P MTCAj) 

until  |TCAj+  |  -  TCAj  I  <  .05  sec 
where  TCAn  =  assigned  TCA  value  (as  defined  above) 
p  =  r^s  r()s  =  range  vector  at  I'CA 

P  -  rNS  f,|s  and  p  =  rNS  -  V„s 

r^s  =  NAVSAl  inertial  position  at  TCA 
rIIS  =  l'os|  satellite  inertial  position  at  ICAj 

Ibis  iteration  also  works  for  type  5  passes.  If  (lie  assigned  ICA  is  negative  (types  I  or(>l.  no 
iteration  takes  place  because  no  real  ICA  exists  for  these  pass  segments.  The  elevation  angle  to  the 
NAVSAl  is  defined  as  the  angle  between  the  range  vector,  p.  and  a  plane  perpendicular  to  the 
vector  from  the  center  ot  the  earth  to  the  host  satellite,  A  positive  elevation  angle  means  the 
NAVSAl  is  above  this  plane. 


Flevation  angle  =  ‘)0  Arecos  ( - — 

\  Ip!  lr,,sl 


The  azimuth  angle  is  defined  as  the  angle  in  the  plane  perpendicular  to  r1)s,  1‘ns-  between  the 
instantaneous  orbit  plane’s  intersection  with  l’||S  and  the  projection  of  the  range  vector,  p,  on  n§ 
measured  clockwise. 


0  =  Arccos 


'(rns  x  rvs 1 "  (rns  x  rns) 

|rtlS  X  rNS  rllS  X  rHS‘ 


().  <0  <  ISO0 


then  the  azimuth  angle  = 


360:  0  if 


0 


rHS  X  rHS 

lrns  x  rns 1 

otherwise 


r\S  >  °- 


In  addition,  at  TCA  an  approximate  longitude  and  latitude  ol  the  subsatellite  point  is  computed 
as  follows. 


Define  AT  =  T(  A  -  t 


VI 


w  =  arctan  (yj|jj/X||S)  -  mod  ^  wAl.360. 


( - co  AT.  360?  \ 


x  =  (  xf,s  +  yfis'1 12  cos  ^ 


y  =  <xns  +  yf,s»l/:si"  * 

Then  the  longitude  and  latitude  are  given  approximately  by: 

A  =  arnanty  xi  0.  <  A  <  3oO. 


0  =  Arc  tun 


'US 

A  A  |  /A 

i(xfi.s +  yfis  * 


00.  K<t>  <60? 


where  tv  (  =  recent  time  ol  vernal  equinox 

w  =  I  artli’s  sidereal  rotation  rate  =  7.2621 15X55  X  10  ^  rad/sec 

rHS  '  host  satellite  inertial  position  at  K  A 

I  he  number  ol  losses  ol  lock  and  reacquisition  tor  each  pass  is  computed  by  counting  the 
number  of  0.  observation  values  and  subtracting  I  lor  the  first  dummy  observation  in  the  pass  and 
subtracting  the  number  of  observations  set  to  0.  because  ol  failure  ol  the  range  difference  absolute 
tolerance  test.  These  losses  of  lock  are  also  summed  over  all  passes.  The  number  of  receiver  switches 
for  each  pass  and  over  all  passes  are  also  computed.  Also  for  each  pass,  the  number  of  z.ero  retrac¬ 
tion  count  observations  is  computed.  In  addition,  the  percent  ol  zero  retraction  count  observations 
for  all  passes  is  determined 


This  pass  through  the  NAVI* AC'  Doppler  data  results  in  a  CELEST-formafted  time-corrected 
observation  tile  with  the  pass  ordering  based  on  the  sorting  procedure.  The  header  record  lor  each 
pass  contains  the  pass  numbers  as  they  appear  in  the  pass  summary  table  (for  use  in  QSOLVE)  and 
tile  pass  type  number  (for  use  in  QSOLVE  and  CELEST ). 


REORDERING  PASSERS  BY  NAVSAT 

The  last  pass  through  the  NAVPAC  Doppler  data  sorts  the  passes  by  NAVSAT,  which  automa¬ 
tically  time  orders  the  data  for  each  satellite  by  TCA  except  for  any  segmented  passes  of  priority 
defeat  data.  In  this  case,  all  the  pass  segments  from  one  receiver  appear  in  sequence  followed  by  the 
same  sequence  for  each  of  the  other  two  receivers.  The  data  are  sorted  by  NAVSAT  because  both 
QSOLVE  and  CELEST  process  all  the  data  from  one  satellite  before  going  to  the  next  satellite. 


PRELIMINARY  DATA  DUALITY  r\  aLUATION  (QSOLVE) 


The  NAVPAC  SS  I'  data  quality  is  evaluated  a  ■  days  after  data  collection  using  the 

QSOLVE  computer  program.  This  program  is  simi.'  u  to  the  CELEST  1  1LTER  section  but  contains 
features  pertinent  to  more  timely  data  analysis.  It  requires  a  CELEST-formatted  perturbed  trajec¬ 
tory'  for  the  host  satellite  and  inertial  trajectories  for  the  NAVSATs  (up  to  6).  In  addition,  the 
program  uses  the  SST  observation  file  created  b/  the  inAVPAC  program.  Only  the  data  correspond¬ 
ing  to  the  overlapping  host  and  NAVSAT  trajectory  spans  are  processed.  The  data  for  each 
NAVSAT  are  processed  a  pass  at  a  time  sequentially.  The  results  for  each  pass  are  summarized  at 
the  end  of  the  run  and  labelled  with  a  pass  number  used  to  relate  the  results  back  to  the  NAVPAC 
program  output. 

QSOLVE  employs  a  least-squares  differential  correction  technique  to  remove  any  signal  left 
in  the  residuals  for  a  pass  (based  on  preliminary  trajectories)  and  estimates  the  variance  of  the  data 
based  on  linearly  adjusted  residuals  after  fit.  The  assumption  is  made  that  the  residuals  can  be 
expressed  as  a  linear  combination  of  the  parameters  to  be  solved  for.  Points  inconsistent  with  the 
rest  of  the  points  are  deleted  and  another  fit  done.  This  iterative  procedure  continues  until  either 
the  maximum  number  of  iterations  is  reached  or  until  two  consecutive  iterations  with  no  points 
deleted  occur.  The  parameters  that  can  be  adjusted  to  remove  any  signal  present  are  frequency  bias, 
frequency  drift,  and  radial,  along-track,  and  cross-track  position  and  velocity  of  the  host  satellite  at 
TCA.  No  parameters  related  to  the  NAVSAT  trajectory  are  present  in  this  adjustment. 

I  or  each  pass  (or  segment  of  a  pass  for  long  passes  that  have  been  split),  the  evaluation 
proceeds  in  the  following  steps. 

I .  Compute  the  TCA  for  the  pass  using  the  iterative  procedure  given  in  the  Data  Preprocessing 
section.  If  the  initial  TCA  is  negative,  no  TCA  iteration  is  possible.  Compute  the  elevation  angle  at 
TCA  given  by  the  formulas  in  the  Data  Preprocessing  section.  Compute  the  rotation  matrix  to  be 


used  to  transform  the  equations  from  the  inertia!  Cartesian  reference  frame  to  the  RAC  (radial, 
along-track.  and  cross-track)  reference  frame  at  TCA  as  follows: 


Let  l!  =  — ~ 

irus| 

rHSXrHS 
W  =  . 

|rl! S  X  rlis| 

V  =  W  X  0 

where  and  r^  are  the  host  satellite’s  inertial  position  and  velocity  vectors  at  TCA 

then  Ktca  =  <U  V  W) 

tx,l 

')rHS  rIIS<  A) 

Also  save  ■C'CTCA)  = - =  6  X  6  state  transition  matrix 

3ellS 

2.  l  or  each  observation  in  the  pass  do  the  following: 

a.  Interpolate  off  of  the  host  satellite  trajectory  to  get  r  n$U0bs.  I  and 
3rHS<  lobs,  I 

^r(tobs  1  =  — T - =  uPPyr  x  fl  portion  of  the  state  transition  matrix. 

‘  aeHS 


Interpolate  off  of  the  NAVSAT  trajectory  to  get  rNS(t„hsjL  To  get  r^s  at  the  time  the  signal  left 
the  NAVSAT  an  iterative  procedure  is  required. 

Let  AtTT  -  lrNS'*o* ”  rHS*to*l  A  *0  ~  ’obsi 
Define  t,  =  t()  -  At^ 

Compute  a  new  transmission  time  estimate 
AtVr  ~  lrNS( * J  rIIS(,o)l  A 

If  |Atjj  -  Aty  j  j  >  10  6  sec,  continue  this  procedure  by  letting  tj+ t  =  tQ  -  Aty|  and 

iterating  until  jAtyy1  1  -  AtyyJ  <  10"*’  sec 


Then  compute  and  save: 


P(,.>bs  >  '  rNSaobs  -  Al T  T  >  _  rUS(W 


f  igure  4  SST  Observability  Geometry 


where 


where 


c.  Assuming  that  the  relevant  quantities  for  the  observation  at  tobs.-,  have  been  saved: 
I )  Obtain  the  computed  range  difference  value: 

D;  =  p*  {obs; '  ~  p(  ,  ’ 


^P^obs.^ 


2)  Compute  the  required  partial  derivatives  of  p  at  t()bs_- 

c  =  speed  of  light  in  km/sec 

re  a> 


9fu 


r„ 


(t 


ohs 


fo=  10 


^P^obs,  * 


=  _  S.  (t  .  -  TCA):  X  86400. 
•  obs. 


fb  in  ppm 
fd  in  ppm /day 


3p(tobs)  3P<W> 

— - 1_  =  - - - ->r(t  obs.)0_I(TCA) 

3rr(TCA)  rlr|IS<tobs>  r  s‘ 

3P(tobs  1  P^obs.* 


^HS^obs,*  lp(t<>b.sj 

3rllS(tohs  1 


^r(tobs»; 


Be 


US 


3)  Form  the  A,  matrix  and  O-  C for  this  observation  as  follows: 

and  O  -  Cj  =  Observed  valuej  -  CRDi 


.  / 

m  3Di 

31),  ^ 

i= 

afb  9|j 

3r  rtTOA ) ^ 

1 

1X8 

''’P^obs,* 

9p(tobs,_,  1 

K 

9»b 

3fb 

3D: 

^P^obs1 

9P^ohSj_|  * 

9f<! 

9f,l 

10 


I 


aDj  "R'W1  ^(tohS). ,» 

dr  rdCAl  r)rr(TCA)  drr(TCA) 

4)  Compute  the  weighted  sum  of  scpiares  of  the  original  O  -  (  s: 

weighted  SOS  = 


3.  If  the  number  of  observations  remaining  is  <  X,  processing  for  this  pass  is  terminated. 
Return  to  step  I.  for  the  next  pass. 

4.  Initialize  the  following  quantities 

Sigma  multiplier  =  amu„  =  I .  Sigma  toleranee  =  o,ol  =  10000.  Ns  =  0 

5.  Do  absolute  and  sigma  tolerance  testing  on  each  point.  (On  the  first  pass  through  this 
computation,  the  original  O  -  Cs  are  tested.  Also  those  observations  with  <7j  =  1000.  are  deleted 
from  further  processing) 

If  |()  -  C|ajj  >  absolute  tolerance,  delete  the  point  from  further  processing 
If  (O  -  C|adj  >  "ml.  "mult.  °i-  ^lc,c  ,hc  P°in( 


£ 


to  -  C) 


>rK, 


where  Oj  is  the  original  observation  sigma  on  the  input  SSI  observation  file 

6.  Again  if  the  number  of  observations  remaining  is  <  X.  processing  for  this  pass  is  terminated. 
Return  to  step  ).  for  the  next  pass. 


7.  Compute  the  estimated  variance  for  the  remaining  pass  data: 


estimated  variance  = 


2  «>-(•£« 

i=  l 

nk  +  na  Ns 


where  =  no.  of  points  remaining 

Na  =  no.  of  parameters  with  a  priori  sigmas 
=  no.  of  parameters  solved  for 


X.  If  zero  points  were  deleted  on  this  iteration  and  also  on  the  previous  iteration,  go  to  step 
lb.  Also  increment  the  iteration  number,  and.  if  greater  than  the  maximum,  go  to  step  lb.  Other¬ 
wise  continue  with  the  next  step. 


9.  Recompute  the  sigma  tolerance 


"t„l.=  f:  + 


N..  +  50./  N 


N„ 


R 


where  N0  =  no.  of  observations  present  after  the  observability  test  was  completed. 

10.  Compute  the  pass  normal  equations  for  the  remaining  observations  using  their  most 
recently  estimated  variances  as  follows: 


BAp  =  l 
=  aT 


where  B  =  A  WA  A  =  [  : 


A> 


AN, 


()  -  c  = 


O-  c, 


O  -  c\ 


and 


l  =  AtW(()  -  O 


org. 


W  =  diag  {- 


'  *ai  ^tnult 1 


i  =  I.: . Nr 


1 1 .  Transform  these  normal  equations  so  that  the  position  and  velocity  adjustments  are  in 
the  RAC  coordinate  system  instead  of  the  TCI  coordinate  system  at  TC'A: 


Define 


wlu'ie 
and 

12.  Add  the  a  priori  parameter  information  into  BR  :  i.e..  add  I  /ct^  to  each  diagonal  for  which 
On  -A  ()  Do  a  singular  solution  of  the  resulting  mutual  equations  using  the  PASMAT  routine  from 
(11  l.SI  (see  Reference  5)  to  get  ApR  and  BR'  ApR  will  have  zero  rows  and  BR'  will  have  zero 
rows  and  columns  corresponding  to  parameters  not  solved  for. 

Compute  Ap  =  R  ApR 

and  B" 1  =  R  BR*  Rr 


BR  APr  =  1  R 

Bk  =  RrBR 
1r=RtI 
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> 


13.  Compute  the  weighteil  sum  of  squares  for  the  points  used  in  the  above  solution 


N 


un¬ 


weighted  SOS  =  23 

i=1  ((,i°,milt>" 


anil  a  linear  adjustment  to  this  quantity  1  r  Ap .  Then,  define  the  predicted  signal  to  noise 


S/N 


Preil. 


weighted  SOS  -  I 

N.<  +  \\  '  N 


Ap 


i/: 


s 


where  Ns  =  no.  of  parameters  solved  for  in  I’ASMAT 
Redefine  omu„  as 

°mult.,new>  =  S/NPred.  X"mul..(old) 

14.  Linearly  adjust  the  O  -  Cs 


<°-C)adj  =(°-Oorg.-  AAP 

and  compute  an  estimate  of  the  variance  before  tolerance  testing 

Nr 

l.stimated  variance  =  1  1  ■  - - — 

nk  +  na  -  Ns 


15.  Return  to  step  5.  to  continue  the  processing. 

16.  Processing  for  the  current  pass  has  now  been  completed.  Information  for  a  summary 
table  is  saved  and  step  I.  is  initiated  for  the  next  pass. 

Various  levels  of  printout  arc  available  during  each  iteration  above.  The  most  useful  is  a  printer 
plot  of  the  adjusted  O  -  Cs  in  which  the  horizontal  spacing  between  points  is  uniform  even  if  the 
observations  are  not  uniformly  spaced.  Any  signal  left  in  the  adjusted  C)  -  Cs  is  usually  apparent  in 
these  plots. 


A  summary  table  is  given  after  all  passes  have  been  processed.  It  contains  the  following  informa¬ 
tion  grouped  by  NAVSAT  and  ordered  by  TCA  for  each  NAVSAT: 

I.  Two  pass  numbers  the  first  is  the  pass  number  assigned  in  QSOLVL  and  the  second  was 
assigned  in  the  NAVPAC  program. 


! 


i 


i 


2.  No.  of  iterations 


3.  No.  of  points  deleted  by  the  SST  observability  test 

4.  No.  of  points  deleted  hy  tolerance  testing 

5.  No.  of  points  remaining  in  final  iteration 
b.  '*<  of  points  deleted  by  tolerance  testing 

7.  Receiver  number(s)  (contained  on  the  SST  observation  file) 

8..  TCA  (day.  seconds)  and  elevation  angle  at  TCA 

9.  Maximum  ()  -  C "adj,  and  square  root  of  the  estimated  variance  based  on  the  final  iteration 
adjusted  O  -  Cs  (labelled  RMS). 

10.  Position  navigations  in  the  RAC  reference  frame  and  their  sigmas  (from  Apr  and  Br1  ) 
(both  are  0.0,  if  the  parameter  was  not  solved  for) 

1  I .  frequency  bias  solution.  fh.  in  ppm 

1  2.  f  requency  drift  solution,  fj.  in  ppm/day 

13.  NSWC  satellite  no 

ORBIT  DETERMINATION  (CELEST) 

GENERAL 

After  the  NAVPAC  SST  data  have  been  preprocessed,  it  is  in  the  proper  form  for  use  in 
CELEST  NSWC’s  orbit  computation  program.  CELEST  employs  a  classical  weighted  least- 
squares  differential  correction  technique  to  fit  satellite  initial  conditions  and  force  and  measure¬ 
ment  model  parameters  to  various  types  of  observations.  The  program  consists  of  four  major 
sections: 

1 .  ORBGf.N  generates  the  reference  trajectories  and  dynamic  partial  derivatives 

2.  EILTER  edits  and  determines  weights  for  the  data  and  forms  pass  matrices 

3.  BSOLVR  expands  and  combines  pass  matrices  to  obtain  solutions  and  computes 
diagnostics 

4.  COVAR  propagates  the  solutions  and  their  covariances  to  produce  the  fitted  trajectory 

CELEST  is  unlike  most  orbit  computation  programs  in  that  it  employs  the  “pass  matrix”  concept. 
Pass  matrices  are  essentially  normal  equations  based  on  data  from  each  pass  separately.  The  FILTER 
forms  the  pass  matrices  as  an  integral  part  of  its  data  editing  and  weight  determination  procedure. 


Il 


In  the  BSOl.VR  section,  the  solution  for  a  given  span  (either  a  long  or  short  are)  is  haseil  on 
pass  matrices  with  I  (  As  in  this  span  only  I  hese  so-called  canonical  pass  matrices  are  updated  to  the 
epoch  of  the  fit  span  and  then  certain  expansions  have  to  he  performed  to  adjust  the  equations  to 
reflect  the  actual  drag  and  thrust  profiles.  Reference  4  contains  a  complete  mathematical  descrip¬ 
tion  of  the  Cl  1.1  SI  program  before  it  was  modified  to  handle  the  NAVI* AC-related  data  processing 
procedures. 

I  lie  purpose  of  the  rest  of  this  report  is  to  describe  the  modifications  and  additions  to 
Cl  1  I  SI  that  were  specifically  designed  to  refine  the  processing  procedures  for  station  tracking 
data  for  the  host  satellite  and  or  the  VAVSA  I  s  (single-satellite  model  and  to  make  possible 
processing  of  the  N  AVI’AC  SS  I  data  ( multisatellite  mode  I.  Modifications  that  affect  the  NAVSA  I 
processing  were  required  because  the  precise  NAVSA  I  trajectories  (in  the  form  of  hybrid  trajec¬ 
tories)  are  required  for  the  N  AVI’ AC  SSI  data  processing  procedures. 


INTEGRATION  (ORBC.EN) 

All  sections  of  (  I  LI  S  I  were  modified  to  accomodate  up  to  eight  satellites  simultaneously 
the  host  satellite  and  up  to  seven  NAVSA  Is.  ORBGI  N  has  been  expanded  to  integrate  up  to  eight 
trajectories  of  the  same  type  (inertial,  perturbed,  or  earth-fixed )  in  one  run  done  sequentially.  If 
the  initial  conditions  for  the  NAVSATs  are  not  given  at  the  same  time,  each  set  can  be  updated  or 
backdated  to  a  general  epoch  defined  as  the  host  satellite  trajectory  epoch.  The  only  restrictions  on 
these  integrations  are  that  all  NAVSA'I  integrations  must  use  the  same  gravity  field  model  and  are 
limited  to  two  drag  segments  and  one  thrust.  The  primary  use  for  this  capability  is  for  simulation 
studies.  In  practice,  all  NAVSA  I  trajectories  required  for  NAVI’ AC  data  preprocessing  and  (  ELI  ST 
processing  are  determined  independently  by  running  Cl  LI  ST  in  the  single-satellite  mode  using 
station  tracking  data. 

Perturbed  trajectories  ( time  histories  of  position  and  partials  of  position  with  respect  to  initial 
conditions  and  force  model  parameters)  are  required  in  Cl  LI. SI  for  all  satellites  for  which  data  are 
to  be  processed.  Hybrid  trajectories  (modified  perturbed  trajectories  defined  under  the  COVAR 
section  below)  are  actually  required  for  the  NAVSA  Is  when  processing  NAVP.AC  SST  data.  The 
epochs  of  the  NAVSA'I  trajectories  must  be  equal  to  or  earlier  than  the  host  satellite  trajectory 
epoch.  In  the  multisatellite  mode,  if  the  solution  has  not  converged  and  reintegration  using  improved 
initial  conditions  from  the  fit  is  required,  only  the  host  satellite  trajectory  is  reintegrated.  Any 
solved-for  corrections  to  the  NAVSA  I  initial  conditions  are  ignored. 

ORBGI  N  is  also  used  to  generate  the  host  satellite  perturbed  trajectory  for  use  in  the 
NAVI’AC  program  "Quick  Look”  route  and  QSOLVL.  Lor  this  the  ARIK  jusq  density  model 
was  incorporated  into  Cl  LI  S  I  .  This  model  is  defined  as  follows 

l  or  altitudes  7b.  ■-  h-  10X.  nmi: 
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where 


P 

P76 

h 

*10.7 
COS  W 


density  in  slugs./ft 2 
density  at  76  rmii  in  slugs  ft2 

satellite’s  geocentric  altitude  above  its  subpoint  in  ntni 
solar  flux  at  10.7-cm  wavelength 


-  -j  [<*«.  +  y m s )  cos  0  +  (yl’s  -  xms)  sin  0  +  zn^j 


where 


r 


d 

I 


satellite  position  vector  = 


cos  Xs 
sin  Xs  cos  I 
sin  A .  sin  I 
0.5  5  rail 

nd  \ 

182.625  / 


1.41  +  .0335  sin 


(|r|-  =  x-  +y-  +z~) 


trd  \ 
182.625/ 


no.  of  days  elapsed  since  3  I  December  1  ‘>57  to  epoch  day  ( I  yr  =  365.25  days) 
orbital  inclination  (in  rad) 


l  or  altitudes  108.  <  h  <  378.  ntni 

p  =  p0(h)  .85  I'  |07  [l.  +  .02375  (e01021'  -  I.1))  (  1.0  + cos*)  ] 
where  p  ,  |u  =  e2.3<):58.5()<)3<  15.7.™  -  .003681.  +  6.363c'  ",I4XI,I 

These  computations  result  in  densities  in  slugs  1C  and  need  to  be  converted  to  kg/km2  by  multipli¬ 
cation  by  5.15375  X  10' 1  Th  is  model  can  be  used  for  determining  the  drag  accelerations  to  be  used 
in  the  integration  of  the  equations  of  motion  only.  The  variational  equations  were  not  modified  to 
accommodate  this  model:  i  e  .  ftp.'ii r  is  still  derived  from  the  standard  (FUST  density  model. 


DATA  EDITING.  WEIGHT  DETI  RMINATION.  AND  PASS  MATRIX  FORMATION  (FILTER) 
General 

The  purpose  of  t he  <11  I  S  I  I  II  I  I  K  section  is  to  edit  and  assign  weights  to  the  data  based 
on  its  noise  variance  and  to  form  pass  matrices  lor  use  in  IISOl  VR.  l  or  the  NAVI* AC  SST  data 
((lass  4.  I  ypc  2  i.  basically  the  same  piocedutes  as  deli  lied  hi  Reference  3  for  ground  station  tracking 
((■I)  data  are  used  with  modifications  made  to  account  for  the  fact  that  the  data  are  collected  by 
an  orbiting  satellite  instead  of  a  ground  station.  These  modifications  are  given  below.  The  FILTER 
section  can  process  ( ■  I  data  only  .  SSI  data  only,  or  both  (done  sequentially  ( TI  then  SST)  and 
generate  the  appropriate  pass  matrices  Modifications  that  affect  the  (i  f  data  processing  procedures 
(for  both  the  NAVSA  I  s  and  the  host  satellite)  were  included  as  part  of  the  NAVI’AC  changes. 

I  hese  are  described  next 


CiT-Related  Modifications 


The  primary  (I  f-related  modification  was  adjustment  of  the  ptocessing  procedures  to  account 
for  the  fact  that  the  transmitting  antenna  is  displaced  front  the  center  of  gravity  for  each  satellite, 
it  is  assumed  that  the  satellite  is  oriented  with  its  negative  vertical  axis  through  its  center  of  gravity 
(c.g. )  pointed  at  the  center  of  the  earth  and  no  rotation  about  this  axis  (yaw  motion)  is  possible. 

I  his  is  an  approximation,  since  the  NAV'SA  I  and  host  satellite  stabilization  systems  keep  this  axis 
oriented  normal  to  the  ellipsoidal  eartli  using  gravity  gradient  and  active  control  techniques,  respec¬ 
tively.  The  antenna  offset  from  the  c.g  must  be  taken  into  account  when  determining  the  computed 
value  for  each  observation  and  its  partial  derivatives.  The  reference  trajectory  used  is  assumed  to  be 
the  trajectory  of  the  center  of  gravity  of  the  satellite. 

for  range  difference  data  tv  pcs  ((  lass  l>.  Types  b.  *).  and  5 ).  the  inertial  position  of  the  trans¬ 
mitting  satellite  at  the  time  of  signal  emission  is  adjusted  as  follows: 


Compute  a  rotation  matrix  RK^(  by 


kkac 
.)  x  1 


/r  (rXf)Xr  rXf\ 

( - r - - r  =  <C  v  w 

\  Irl  !(f  X  r>  X  n  Ir  X  r|  / 


where  r.r  =  rc  ,  r  =  satellite's  position  and  velocity  vectors  at  time  of  signal  emission 
This  defines  the  assumed  body-fixed  axes  orientation  in  inertial  space. 

fhen  rant.  =  rcg.  +  KRC  AArant. 

where  Arant  =  antl-‘nna  offsets  with  respect  to  the  c.g.  in  the  body-lixed  radial,  along-track.  and 
cross-track  (RAC)  directions 

In  addition,  the  inertial  velocity  of  the  transmitting  satellite  for  Doppler  frequency  data  (Class  7. 
Type  ".  NAV’SA  Is  only  i  is  adiusted  as  follows 

fant  =  rc.g.  +  w  X  HK  V(  Aranl. 


1  ,  .  ',ri 

where  w  =  ~~  tV  •  ri\V  +  ,  (W  -ril  =  KA(  trame  angular  velocity  vector 

Irl  ir  x  r, 

r  =  satellite  acceleration  vector 

Another  modification  involves  the  use  of  temperature  and  pressure  values  as  functions  of 
station  height  as  input  to  the  modified  llopfield  tropospheric  refraction  correction  model.  If  the 
pressure  and  relative  humidity  on  the  observation  file  header  for  a  given  pass  are  both  zero,  the 
following  formulas  are  used  to  define  the  weather  parameters: 

Temperature  I  <  I  =  15.  (>.5  It 

Pressure  (millibars)  =  11)13. 25e 
Relative  humidity  =  XO : 


where  h  =  station  height  above  the  reference  ellipsoid  in  km 


SST -Related  Mollifications 


The  general  How  of  the  N  AVI’AC  SSI  data  processing  procedures  in  the  I  ILTI  R  section  is 
the  same  as  for  station  tracking  data.  Data  editing  and  weight  determination  are  an  integral  part  of 
the  pass  matrix  generation.  Since  the  NAVI' AC  SSI  observation  file  has  the  data  time  ordered  by 
N.AVSA  all  data  for  one  \AVSAT  is  processed  consecutively,  bach  SST  pass  consists  of  a  header 
with  identifying  information  anil  observations  in  the  form  Jtime.  observation,  sigmaj. 

l'he  IT  A  iteration  for  each  pass  is  identical  to  that  described  in  the  Data  Preprocessing  section 
of  this  report  as  are  the  computations  of  the  zenith  angle  (‘>0 -elevation  angle)  and  azimuth  angle  at 
l(  A.  Another  computation  to  be  done  for  t(  A  is  the  formation  of  a  rotation  matrix. 
required  to  convert  from  the  inertial  reference  frame  to  the  navigation  reference  frame. 

/  r  (r  x  r)  *  r  r  x  r 

r<  A  \ trj  Hr  x  f i  x  r|  ir  X  f| 

t  x  .1 

where  r.  r*  =  Pj JS»  IT  At.  i((S(  1<  At 
Also  the  two  (1  X  ti  state  transition  matrices 

/*HS«'r< 'A>\  |  /vNS<l(  Al\ 

\  H  A )  /  \  C  I  (  A )  / 

need  to  be  saved  for  use  in  transforming  the  pass  normal  equations. 

for  each  observation  in  the  pass,  the  following  items  are  peculiar  to  the  NAVI’AC  SSI  data 

1 )  No  station  coordinates  or  tropospheric  refraction  correction  computations  are  applicable 
to  the  SS  I  data  This  also  implies  that  no  AIM'D  (earth-fixed  to  inertial  transformation  I  matrices 
are  required  since  all  computations  involve  trajectories  that  are  given  in  the  inertial  reference  frame 

2i  Observations  with  <>  =  1000  are  deleted  and  counted  as  deweighted  points  even  though  a 
non-zero  observation  value  may  be  present  I  liese  correspond  to  zero  refraction  count  observations. 

T  >  I  he  SSI  observability  test  as  described  in  the  Preliminary  Data  Quality  I  valuation  section 
of  this  report  replaces  the  zenith  angle  tolerance  test 

4)  l'he  time  transmission  correction.  At  j  j  is  applied  at  the  transmitter  (N.AVSA  l  l  as 
described  in  the  same  section  mentioned  in  Ti  above. 

5 )  If  the  host  satellite  data  reduction  is  being  done  using  the  WT1S-72  station  coordinates,  then 
an  adjustment  to  the  N.AVSA  I  \  inertial  position  is  required  after  the  time  transmission  correction 
is  applied.  I  his  is  necessary  because  the  N  AVN  A  I  trajectories  are  computed  using  NWI  -d/  station 


coordinates  that  differ  from  the  W(iS-72  coordinates  hy  0.26"  of  arc  in  longitude.  The  adjustment  is 
done  as  follows: 


cos  0  -  sin  0  (I 

fvS  =  I  sin  0  cos  0  0 

\  0  0 


rNS  =  RNWI  -UZ  ■  W(iS-72rNS 


where  0  -  7.222  X  10  5  deg  =  0.26” 

These  adjusted  NAVSAT  coordinates  are  then  used  in  all  further  observation  processing. 

6)  I  he  antenna  offsets  for  each  satellite  are  accounted  for  in  obtaining  the  computed  observa¬ 
tion  value  and  partial  derivatives  of  the  data  with  respect  to  satellite  coordinates,  f  or  the  host 
satellite,  define  the  rotation  matrix  by 


rhs 


'US 


,rHS  X  fils’  X  r! IS  rIIS  X  rl(S 


|rlis)  j*r||S  X  fils'  X  fiisl  |rHS  X  ’rus| 


where  rn^.  r^  are  host  satellite  position  and  velocity  from  the  trajectory  (i.e..  for  the  c.g.)  at  time 
of  signal  reception  (t()j,s ) 


then 


filS  "  rl IS  +  KHS  ArIIS  nt 

jni  ant. 


where  r(,s  =  inertial  position  of  the  IIS  receiving  antenna  at  t0^s 


•nil 


Arl)S  -  antenna  offsets  wrt  the  c.g.  in  body-fixed  radial,  along-track  .  and  cross-track 
•,n'  directions 

f  or  the  NAVSAT,  define  the  Riy$  rotation  matrix  by  the  same  formula  with  r^g  and  r^s  replacing 
rns  and  f||S  where  these  are  defined  at  the  time  of  signal  emission  ( tt)y>s  -  Atyy  1 


NS, 


rNS  .  +  RNS  ArNSanI 

auj.  ant . 


where  rNS  =  inertial  position  of  the  NS  transmitting  antenna  at  tohs  -  AtT( 


Ar^  =  antenna  offsets  wrt  the  c.g.  in  body-fixed  radial,  along-track.  and  cross-track 
•‘r’1  directions. 


Ihe  radial  NAVSAT  antenna  offset  can  also  be  used  to  adjust  the  NAVSAT's  position  to  account 
for  the  possible  difference  between  the  central  force  term  used  in  the  NAVSAf  orbit  fit  and 
used  in  the  host  satellite  orbit  fit. 


I  bis  offset  is  given  by 


a\S  An 
^rNS  ~  V 

T.uiial  3  ^\S 

where  a^s  =  semimajor  axis  of  NAVSAT  orbit  in  km 
An  =  Pj(s  n\^  |!l  •s'11'  sv‘s'" 


7  )  I'he  observation  equation  for  the  NAVI’ AC  SSI  range  difference  (  Rl>)  data  is  given  as: 


Rl)  = 


p  +  -  fh(t-  TCA I 


obs 


J  'obs  -  -At 

where  p  =  |rNS  ( t  At  j  r  )  -  r(jS  (t)|  in  km 


f(>  =  I0(1  (gives  fb  units  of  ppm  f 
l‘b  =  frequency  bias  (nominally  zero) 

•it  =  range  difference  time  interval  in  seconds  (t()hs  -  At  is  the  observation  time  of  the 
previous  observation) 


The  required  partial  derivatives  for  p  are  given  as. 
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IIS  orbit  or  force  model  parameter 


NS  orbit  or  force  model  parameter 


where 


arIIS 

-  and 

dpk 


,)rNS 

dqk 


are  obtained  from  the  appropriate  perturbed  trajectory  at  the  specified  times.  The  MS  and  NS  epochs 
may  be  different  however. 


Then  the  partial  derivatives  for  RD  are  given  as: 
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A  row  of  the  A-matrix  (partials  of  data  with  respect  to  parameters)  has  the  parameters  (columns) 
ordered  as  follows. 


(p  q  f.)  maximum  of  23  parameters 


where 


aT  qT 
,T  =  /.T 


P  =  <eHS( 


WlQ  ^  « 


qT  =  (e^s  C  ^  T^s  ) 

and  drag,  thrust,  and  radiation  pressure  parameters  are  optional. 

After  all  the  observations  for  a  pass  have  been  processed  as  described  above,  the  A  and  ()  -  ( 
matrices  for  all  observations  are  converted  to  the  normal  equation  matrices  B  and  I  as  follows: 


B  =  A 1  W  A 


=  AT  W<<)  -  (  ) 


where 


W  =  ilia  e 
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'obs 


\"ohs  / 

noise  variance  on  the  i'1'  observation 


I  or  the  adjustment  procedure  (required  for  data  editing  and  weight  determination'  a  subset  of  the 
above  normal  equations  is  used  and  is  represented  by 


B  =1 


'h.b 

»h.l 

»b.:\ 
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b.l 

»l.l 
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and  1  „ 

1 

’h,2 

B,  s  / 

13  x  13 

where  subscript  b  refers  to  frequency  bias  ( 1  I 

subscript  I  refers  to  host  satellite  orbital  elements  <(>) 
subscript  2  refers  to  NAVSA  I  orbital  elements  tot 

30 


I'lie  matrix  required  to  transform  these  normal  equations  into  equations  for  inertial  position  and 
velocity  at  TCA  is  given  as 

/ 1  O 

i^( J'C'A)  =  I  O  i^jtK  A) 

13  X  13  \()  O 

where 


The  actual  transformation  is  given  hy 

B(  TCA  i  =  v  "T<  If  A)  B()  0  'lICAl  and  l(H  At  =  v  I C  A)l(> 

To  transform  this  system  of  equations  into  equations  for  IIS  navigation  reference  frame  positions 
and  velocities  at  TCA,  the  following  matrix  and  transformation  are  required: 


/  Vhs*  ,(  A»\ 
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ktca 

O  kT(A 

1 3  X  13 

then  B'tTCAi  =  RrH(T(  A)R  and  I  VlCA)  =  R1  1  ( ICA) 

I  he  a  priori  statistics  on  the  frequency  bias  and  the  six  NAVSAT  parameters  are  included  hy 
adding  I  Up  terms  into  the  appropriate  diagonal  elements  of  B'.  Then  the  equations 

B’(TCA)  Apt  TCA )  =  I  'l  if  Al 

are  solved  using  the  I’ASMA  I  routine  modified  to  accommodate  six  more  parameters.  These 
additional  parameters  are  treated  identically  to  the  six  host  satellite  parameters.  Ihe  Apt  TCA) 
is  then  transformed  hack  to  bias  and  orbital  element  corrections  at  epoch  as  follows: 

Ap()  =  Cc’n'CA)  R  Apt  1(  A) 

Tile  adjusted  O  -  Cs  are  then  given  by 

(°-C)ad|  =tO-(l,)rg  -A„Ap0 
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I  or  the  NAVI* AC  SS  I  data,  the  predicted  signal  to  noise  for  a  pass  is  given  as 
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i 


S  N 


Variance  -  I  Ap„ 

pK'a  =  ^  nk  +  na  Ns 


i/: 


where  \p  =  no.  of  points  remaining 

\  A  =  no.  of  parameters  with  a  priori  sigmas  =  7 
Ns;  =  no.  of  parameters  solved  for  in  PASMA  I 

As  was  mentioned  earlier.  QSOLVT  and  the  (  I  I  I  SI  11111  R  section  are  similar  in  the  way 
each  determines  the  noise  variance  on  the  data  for  each  pass.  The  main  differences  are  listed  below: 

1  )  Tll.TLR  can  also  adjust  NAVSA I  position  and  velocity  components. 

2)  QSOLVT  can  solve  for  an  additional  bias  parameter  frequency  drift. 

3  i  Tacit  iteration  in  l  I  LI  I  R  retests  points  that  were  previously  tagged.  In  QSOLVT.  once  a 
point  is  deleted  it  is  never  reconsidered  in  later  iterations. 

4i  The  sigma  tolerance  can  be  fixed  by  input  in  TILTI  R  but  not  in  QSOLVT. 

5  i  I  lie  estimate  of  data  quality  in  (III  ST  is  called  the  “filtered  Noise"  and  is  defined  as 


1/2 


1  (0,0 


r  mull . 


filtered  noise  =  I 


QSOl  VI  uses  the  adjusted  ()  Cs  to  estimate  the  standard  deviation  of  the  noise  defined  as 


1/2 
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SOLUTIONS  AND  DIAGNOSTICS  I BSOLVR ) 
General 


BSOI  VK.  the  solution  and  diagnostic  section  of  Cl  I  ISI ,  has  been  modified  to  process  the 
N  AVI’AC  SS  I  data  separately  or  simultaneously  with  station  tracking  data  for  the  host  satellite  The 
program  can  accommodate  parameters  for  up  to  seven  NAVSA  Is  along  with  the  host  satellite 
parameters.  The  total  number  of  arc  parameters  is  limited  to  K7  w  ith  a  maximum  of  3‘*  host  satellite 
arc  parameters  allowed 

If  station  data  pass  matrices  are  present,  they  are  time  updated,  expanded,  bias  eliminated,  and 
summed  first  If  \AVT’A(  SSI  data  pass  matrices  are  also  present,  they  are  processed  similarly  for 

,?2 


\ 


each  NAVSAT,  ami  then  combined  with  the  (i  I  pass  matrices  to  obtain  tile  solution  for  all  arc 
parameters  simultaneously.  1  he  epoch  for  all  arc  parameters  is  defined  as  the  host  satellite  trajectory 
epoch  for  long-arc  fits  and  the  program  determined  epoch  for  short-arc  fits,  l  or  this  reason,  the 
NAVSAT  trajectory  epochs  must  be  the  same  or  earlier  than  the  host  satellite  trajectory  epoch. 


Two  modifications  to  BSOl.VR  were  implemented  to  improve  the  long-arc  fit  diagnostic 
capabilities  of  the  program  Ihe  first  was  an  option  to  compute  navigations  before  fit  without  doing 
cross-pass  filtering.  I  he  second  was  an  option  to  cross-pass  filter  based  on  ratios  of  navigations  to 
their  estimated  uncertainties  instead  of  navigations.  This  ratio  is  called  a  “signal  to  noise"  and  is 
given  by 

-^ri 

SNlAr^  =  -  —  i  index  signifies  a  particular  navigation  component 

where  Ar^  =  navigation  value 

flAr  =  formal  navigation  standard  deviation  from  the  covariance  matrix  of  the  solution 
%rhit=  input  a  priori  component  uncertainty 

In  addition,  all  RWS  navigations  are  then  computed  as  follows. 
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j  index  signifies  a  particular  pass 


The  same  a  priori  orbit  sigmas  are  used  for  station  and  SST  cross-pass  filtering  even  though  their 
navigation  reference  frames  have  only  the  host  satellite  along-track  direction  in  common. 


I  lie  SS  I-related  modifications  are  given  next  followed  by  an  alternate  pass-matrix  expansion 
technique  for  the  host  satellite  parameters  called  the  multiple-velocity  expansion  method. 


SST-Relaled  Modifications 


Adjustment  of  the  SS  I  pass  matrices  to  reference  the  fit  epoch.  ts.  and  expansion  to  accom¬ 
modate  the  drag  and  thrust  segmentation  are  done  as  follows. 

l  et  B(t„ I  and  l  it,,)  represent  the  canonical  pass  matrix  and  its  right-hand  side  for  a  particular 
pass.  Updating  these  normal  equations  to  he  for  orbital  elements  at  ts  results  in  the  following 
matrices. 


lllsT'c'  _ 


‘K'ns<,sl 
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dr\S  rNS<’s) 
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(■>  X  (>  matrices  obtained  from  the  PARDl'L  subroutine 

olCl  I  I  ST 


lj  ami  li  are  identity  matrices  of  different  dimensions 

l  lte  expansion  to  accommodate  the  drag  and  thrust  profiles  and  the  radiation  pressure  update  for 
the  host  satellite  is  done  b\  treating  all  the  NAVSAT  parameters  and  the  frequency  bias  parameter 
as  the  bias  parameters  are  treated  presently  in  the  station  pass-matrix  expansion  procedure  (see 
Reference  3).  After  this  is  completed,  the  same  procedure  is  used  to  expand  the  NAVSAT  dynamic 
parameter  sections  only  the  expanded  host  satellite  parameters  and  frequency  bias  are  now  treated 
as  biases.  The  formulas  for  each  of  the  expansions  above  are  a  function  of  the  SST  pass  Tf  A  time 
relative  to  the  drag  and  thrust  segment  times.  Tor  the  NAVSATs.  a  maximum  of  two  drag  segments 
and  one  thrust  is  allowed. 

The  expanded  pass  matrix  after  the  a  priori  frequency  bias  information  is  included  has  the 
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where  B_  K  =  Bit.)  +~ 

n'b  s  b.b  or 
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The  next  step  is  to  formally  eliminate  the  frequency  bias  parameter  from  the  corresponding  equa¬ 
tions.  This  results  in  the  following  matrices: 

Belim.  =  B<u>  '  Bo.b  Bb.b  Bo.b  ami  '  elim.  =  '  <>  "  Bo.h  Bb.b  '  b 

These  eliminated  matrices  (containing  arc  parameters  only)  are  then  summed  for  all  passes  for  a 
given  NAVSAT  within  the  fit  span  except  for  those  that  have  been  tagged.  This  results  in  up  to  N 
( no.  of  NAVSATs)  sets  of  normal  equations  given  by: 
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If  station  tracking  data  are  present,  it  would  already  be  reduced  to  a  system  of  arc  normal  equa¬ 
tions  for  the  host  satellite  parameters  given  by 


B'HS  -  1  MS 


I  host-  matrices  arc  then  combiner)  as  follows  to  form  the  are  normal  equations 


15  Ap  =  F 


Next,  the  a  priori  information  on  each  host  satellite  and  NAVSAT  arc  parameter  is  introduced  into 
the  normal  equations  by  adding  terms  of  the  form  1  lo-  to  each  diagonal  element.  Different  r  allies  of 
Up  for  the  same  parameter  but  for  different  NAVSAT s  can  be  used.  15  is  then  inverted  and  the  equa¬ 
tions  solved  to  get  the  are  parameter  solution  Ap  =15"*  E  and  the  covariance  matrix  B  1 . 

Signal-lo-noise  values  (old  and  bias-red  need  I  are  computed  separately  for  (if  and  SST  data  and 
also  combined.  The  formulas  for  the  SS  I  signal-to-noise  quantities  are  identical  to  the  those  used  for 
the  ( •  I  data  The  combined  signal-to-noise  values  are  given  as  follows: 
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where 


S/N(;T<old)  =  signal  to  noise  before  fit  for  (IT  data 
S'Nsst (old )  =  signal  to  noise  before  fit  for  SST  data 
N|  =  no.  of  observations  in  i*1'  (IT  pass 

Mj  =  no.  of  observations  in  i1*1  SST  pass 

=  total  no.  of  (IT  passes 
N«jt  =  t°tal  no.  of  SST  passes 

S/N(y(adj.)  =  bias-reduced  signal  to  noise  for  (IT  data 
S/NSsT<ad'  )=  bias-reduced  signal  to  noise  for  SST  data 

The  predicted  signal  to  noise  for  the  arc  is  then  given  by 

S/N( predicted)  =  ^S/N* ( bias-reduced )  - 
and  the  percent  change  in  variance  is  given  by 

100  ETAp 

Percent  change  in  variance  =  — - 

N(;t  nssi  \ 

Z  Nj  +  Z  Mjl  S/N- ( bias-reduced) 

i=  1  i=  1  / 

Navigations  are  computed  for  each  pass  based  on  the  expanded  pass  matrices.  For  each 
station  pass  matrix,  the  Spn;,  solution  augmented  with  zeroes  for  the  pass  bias  parameters  is  used 
to  linearly  adjust  the  right-hand  side  and  then  the  present  form  3  navigation  procedure  is  followed. 
For  each  NAVPAC  SS  I  pass  matn.x.  the  navigation  procedure  is  as  follows: 

I )  The  sections  of  the  expanded  pass  matrix  corresponding  to  the  NAVSAT  orbital  elements 
and  the  frequency  bias  before  the  bias  elimination  was  performed  are  required  for  the  navigation 
procedure. 


7X7  7X1 


where  '  indicates  only  those  rows  and  columns  corresponding  to  orbital  elements  for  the  NAVSAT 
are  present. 


2)  The  rijzli (-hand  side  is  then  linearly  adjusted  t<>  account  lor  the  solution  lot  the  host 
satellite  parameters  only 


2)  ami  I  p  are  then  transformer!  to  relerenee  NAVSAI  position  ami  velocity  in  the  IIS 
RAC  reference  frame  at  It  A 

HgtTCAl  R'v  '(If  AtlT^  'iK  AiR 

I  ”l  ICAi  =  r'  u  1 1  l<  All  " 

where 

;,r\S  fNs"(  A> 

C/(T(At=  ■'  — 

<)eNS(ts. 

/rt(Ao  o\ 

R=  °  K,(a  O 

\o  <>  >/ 

Rj{A  =  RAC  transformation  at  TCA  for  the  pass  do  final  in  the  I  ll.Tl'.R  section 

4i  The  resulting  system  ot  navigation  equations  BgtTCAlAN  =  1  j't'K  At  is  then  solved  using 
the  P ASM  AT  subroutine  that  iletermines  which  components  (may  include  velocity  components) 
can  be  solved  for.  PASMAT  returns  /ero  values  for  a  component  and  its  sigma  if  it  is  not  solved  lor 

After  the  navigations  are  performed  for  all  SST  passes,  the  RMS  and  RWS  navigations  for  each 
NAVSAT  and  combined  are  computed.  The  RMS  navigation  computations  account  for  the  fact 
that  not  all  navigation  components  are  present  for  each  SST  pass. 

The  frequency  bias  solution  for  each  pass  based  on  the  arc  solution  for  the  host  satellite  and 
NAVSAT  parameters  is  given  by 

^'h  =  Bhh(l  h  ‘  'VlS.h  AP|IS  '  '*NS.b  ^PnS* 

Hie  cross-pass  filtering  has  been  extended  to  include  the  satellite-to-satellite  data.  Separate 
first-order  polynomial  fits  to  the  SST  navigations  from  all  NAVSAI  S  are  done,  in  addition  to  the 
fits  for  the  station  navigations.  The  procedure  terminates  only  when  the  same  station  passes  and 
the  same  SST  passes  are  tagged  on  successive  iterations  or  the  maximum  number  of  iterations  is 
reached.  I  or  the  three  SST  navigation  polynomial  fits,  only  those  navigation  components  that 
are  not  zero  ti  e.,  only  solved  for  components)  are  included  in  the  fits.  This  means  that  the  number 
ol  passes  for  each  polynomial  fit  will  be  different  for  each  component  and  usually  less  than  the  total 


number  of  untagged  SS  I  passes.  The  navigation  components  determined  lor  a  particular  SST  pass 
do  not  change  from  iteration  to  iteration.  Both  “total”  anil  “solution"  RMS  and  RWS  navigation 
values  are  computed  for  each  iteration  “Total”  navigation  values  are  based  on  all  passes  that  arc  not 
deleted  by  input  and  are  (lie  same  for  all  iterations.  "Solution"  navigation  values  are  based  only  on 
the  untagged  passes  used  in  the  solution  lor  each  iteration.  All  SST  pass  tags  are  then  honored  in  any 
subsequent  long-arc  or  short-arc  fit 


Multiple-Velocity  Pass-Matrix  Expansion 

The  multiple-velocity  pass-matrix  expansion  allows  fora  simultaneous  solution  for  epoch 
orbital  position  and  velocity  and  additional  small  velocity  increments  at  equally  spaced  times  along 
the  host  satellite  trajectory.  The  velocity  increments  are  solved  for  in  the  RAC  reference  frame  at 
each  time.  This  expansion  can  be  used  for  long  arcs  as  well  as  short  arcs.  There  is  no  reintegration 
possible  after  a  multiple-velocity  fit  and  only  propagated  trajectories  after  fit  can  be  computed. 

This  means  that  a  converged  orbit  using  standard  techniques  needs  to  be  obtained  before  this  solution 
is  performed. 

The  arc  is  partitioned  with  N  velocity  increment  time  points -Jt  j .  ts . t*c|  as  shown  below 

t  =  epoch  of  the  host  satellite  trajectory 
ts  =  epoch  of  the  fit  span 
tt  =  end  of  the  fit  span 

| - // - 1 - 1 - 1 - // - 1 - 1 - // - 1 

t0  ls  ‘1  *2 


The  multiple-velocity  pass-matrix  expansion  js  based  on  the  following  relationship: 

dDttl  _  /  r)l)(  t )  aixu\  /  ,  ( 
afitj >  \<jr(ts)  aft ts >  /  v ’s  1 

w  here  t  >  tj  >  ts  >  t() 

l)(t>  =  observation  at  time  t 


(  )•  denote  the  last  three  columns  ol  a  (>  X  o  matrix 


and  v,  (tjl  =  vWtjl  <tsl 

yirl  =  (>  X  h  state  transition  matrix  at  r  = 


dr  fir) 
dr  ft t„> 
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Also  define  for  i  =  1.2 . N 


i  =  1.2 . N 


Rj  =  (IJj  Vj  Wj)  =  rotation  matrix  to  go  from  RAC  frame  at  !j  to  the  inertial  frame 

3  X  3 
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r<  t( ».  r ( t ! »  =  inertial  position  ami  velocity  of  the  host  satellite  at  t| 

Only  observations  that  occur  at  or  alter  time  tj  are  used  to  determine  the  required  velocity 
increment  at  t,  I  he  assumption  is  made  that  it  lor  a  pass  1C  A  >  tj.  all  the  observations  m  the 
pass  are  to  be  used  to  solve  lot  Av( .  the  vcloeit  y  increment  at  t, .  (  his  allows  an  expansion  of  t he 
|iass  matrix  m  BSOI  V  K  instead  ol  including  the  velocity  parameters  m  the  formation  of  the 
\ -matrix  lor  each  observation  in  I  1 1  I  I  R 


I  he  start inyr  point  lor  this  expansion  is  the  completion  ol  the  standard  time  update  and 
expansion  currently  done  in  BSOI  VK  except  lor  one  small  change  the  updated  normal  equations 
involve  position  and  velocity  at  ts  instead  of  orbital  elements,  l-.acli  updated  set  ol  pass  normal 
equations  ( for  (i  I  or  SS  I  data)  lakes  the  form:  ^ 


where  x  indicates  the  host  satellite  position  and  velocity  parameters  at  ts 
y  indicates  all  other  parameters 


It  t,  <  ICA  ''  tjt  |  .  then  the  velocity  parameter  expanded  normal  equations  take  the  following 
form 
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where  VhR-  =  (())  lor  j  =  i  +  I ,  i  +  2 .  N 

6  X  d 
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Alter  this  expansion  is  completed  tor  a  particular  pass,  the  biases  are  formally  eliminated.  All  of  the 
resulting  expanded  pass  matrices  for  the  fit  span  are  then  summed  and  combined  (if  SST  data  are 
present)  to  form  the  arc  normal  equations.  A  priori  sigma  information  on  each  velocity  increment  is 
included  before  these  equations  are  solved  by  adding  in  terms  of  the  form  l/tr^to  each 
diagonal  element  Separate  a  priori  sigmas  for  radial,  along-track.  and  cross-track  velocity  increments 
are  included.  All  further  computations  m  BSOLVR  account  for  these  additional  host  satellite  orbit 
parameters. 


PROPAGATION  (COVAR) 

This  section  ol  Cl  MSI.  which  propagates  the  corrections  solved  for  in  the  BSOLVR  section 
and  applies  them  to  the  reference  trajectory  to  get  the  fitted  trajectory,  still  operates  in  a  single¬ 
satellite  mode.  Therefore  only  the  portion  of  Sp  and  B“ 1  from  BSOLVR  corresponding  to  the  host 
satellite  are  required.  Several  modifications  to  this  section  have  been  made  to  accommodate  the 
NAVPAC  data  processing  procedures  and  are  given  below. 

I )  The  covariant  ephemeris  was  changed  to  include  the  earth-fixed  velocity  vector  and  a  full 
6  X  b  earth-fixed  covariance  matrix  (position  and  velocity)  at  each  trajectory  timeline. 

rf.;F  =  (ABCD)  j^r,  -  (CI))TW(CD)r,  J 

C,.F  =  H  C,  Ht 

6X6 

V 

where  rF  rf,  C'|  =  inertial  position,  velocity,  and  covariance  matrix 

I)  =  precession  transformation  matrix  , 

C  =  nutation  transformation  matrix  (  all  ^  X  3 

B  =  earth  rotation  transformation  matrix  f 

A  =  polar  motion  transformation  matrix  j 


/0  -2  0\ 

W  =  (  2  0  o  j 

\  0  0  0  / 

2  =  . 7 2 *■> 2  I  I  5H55  X  It)'4  rad/sec  =  harth’s  mean  sidereal  rotation  rate 

(ABCI)  O 
- ABWC'I)  ABCD 

2)  An  option  to  apply  the  rotation  matrix.  Rnwl~9Z-»WGS-72-  defined  'n  •he'  FILTER 
section  to  both  the  propagated  inertial  positions  and  velocities  at  each  trajectory  timeline  has  been 
included  in  COVAR.  This  is  required  if  the  NAVPAC  orbit  fits  are  carried  out  using  NWL-dY  station 
coordinates.  In  this  ease,  the  rotation  would  not  have  been  applied  to  the  NAVSAT  trajectory  when 
filtering  the  SST  data. 


i 
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3)  The  option  to  output  the  trajectory  of  a  reference  location  on  the  satellite  instead  of  the 
c.g.  has  also  been  included  in  COVAR.  If  irre|-  represents  the  coordinates  of  this  reference  loca¬ 
tion  in  a  body-fixed  reference  system,  then  the  transformation  matrix  to  go  from  this  system  to  the 
inertial  reference  system  is  defined  by 


/  r  irXflXr  rXr  \ 
\  |r|  |<r  X  b  x  r|  |r  X  r|  / 


=  ( U  V  Wt 


where 


r  =  propagated  inertial  position  and  velocity  of  the  host  satellite  at  the 
trajectory  timeline 


I This  RAC  reference  frame  is  the  same  as  the  one  used  to  define  antenna  offsets  in  the  I'lLTliR 
section. ) 


Then 


rref.  fc.g.  +  ^ ref.  ret. 


is  the  inertial  position  of  the  reference  location. 

It,  in  addition,  the  inertial  velocity  is  required  in  COVAR.  an  adjustment  to  account  for  the 
reference  location  offset  from  the  e  g.  is  also  required. 


'ref. 


=  r  +  u  X  R  fir  r 
c.g.  "'ret.  ‘ret. 


where 

1  ,  Irl 

oj  =  —  tV  •  rlW  +  - r  <W*r)U  =  RAC  frame  angular  velocity  vector 

Irl  |r  X  r| 

r  =  satellite  acceleration  vector 

The  usual  reference  location  would  be  the  electrical  phase  center  of  the  transmitting  antenna 
l used  for  the  NAVSATs).  However,  the  radial  reference  location  component  can  also  be  used 
to  account  for  differences  between  the  g  value  used  in  the  fit  and  a  desired  value.  For  instance 
if  the  NAVI’AC  orbit  fits  were  done  in  the  NWL-dZ  reference  system  (that  used  for  the  NAVSATS) 
and  the  output  trajectory  was  required  in  the  W(;S-72  system  (that  used  for  the  host  satellite)  the 
radial  of  fset  would  be  defined  as  follows: 

Ar  -  aiis  Am 
Ar"s™t,>  ‘  T 

where  a|]s  =  semimajor  axis  of  the  host  satellite  orbit  in  km 
Am  =  is  -  f'NS  i”  kin-* /sec- 

4)  The  NAVPAC  SST  data  fits  require  a  special  form  of  NAVSAT  trajectory  called  a  hybrid 
trajectory.  A  hybrid  trajectory  is  identical  to  the  original  perturbed  trajectory  except  the  positions 
have  been  updated  to  reflect  the  results  of  an  orbit  fit.  T  he  partial  derivatives  remain  unchanged. 
COVAR  was  modified  to  create  a  hybrid  trajectory  on  option. 
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5)  To  get  a  propagated  trajectory  ;iml  covariance  matrix  when  the  multiple-velocity  expansion 
is  used  in  BSOl.VR.  additional  changes  to  TOVAR  were  required.  I  igure  5  gives  a  representation  of 
the  actual  corrections  solved  for  in  the  multiple-velocity  method,  l-.ach  Avj  is  actually  an  additional 
correction  to  the  inertial  velocity  correction  at  tj  beyond  that  propagated  from  all  previous  2v;’s. 

i  =  I . i  -  I  and  Arj^  and  Arj|c; .  This  results  in  a  trajectory  that  is  continuous  in  position  but 

discontinuous  in  velocity  at  times  tj.  i  *  1.2 . N.  The  state  transition  matrix  required  to  propagate 

the  correction  vector  Apji^t’^)  and  its  covariance  matrix  I5j^ils)  to  position  and  velocity  corrections 
ami  their  covariance  matrix  at  lime  1  is  given  In 

'R,  to  -  (  yi,  <tl  v,  mV,  R, 

6  X  (6  +  dN  +  S  .  l  (<  x  (1  6  X  X 

u 

VjRj  =  (O)  il'l'-t, 

fi  X  .i 

where  the  V.  and  R,  matrices  are  defined  in  the  BSOL.VR  section 


where  d  represents  all  other  dynamic  parameters  (e.g.  drag,  thrust, 
or  radiation  pressure) 

Nj  =  no.  of  these  other  dynamic  parameters 
/  Aril)  \ 

Then  .  ^(UlSPusIt,) 

\  Art  t )  / 
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figure  5.  Multiple-Velocity  Corrections 
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